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Interim  Report  on  Grant  75-2812 
"Numerical  Studies  in  Computer-Aided  Design" 

D.  A.  Calahan,  Principal  Investigator 
University  of  Michigan 
Grant  Period:  4/1/78-3/31/79 

I.  Research  Objectives  in  Progress 

Two  study  objectives  were  supported  by  this  grant. 

A .  Study  of  vectorizability  of  aerodynamic  fluid  flow  codes 

(AFFDL ) .  An  explicit  3-D  aerodynamic  simulation  code  (ex¬ 
ercised  in  2-D )  was  vectorized  in  Fortran  and  coded  in  CRAY-1 
assembly  language.  The  latter  was  greatly  assisted  by  use 
of  a  CRAY-1  simulator  and  a  cross-assembler  developed  under 
grant  sponsorship  [4] {5].  Major  subalgorithms  from  these 
codes  were  benchmarked  and  reported  in  [3].  The  complete 
code  has  also  been  vectorized  in  Fortran  and  assembly  coded 
and  is  planned  for  tinal  benchmarking  this  spring.  Speeaups 
on  the  order  of  100  are  indicated  viz-a-viz  the  C DC  6600. 

E .  Study  of  the  development  of  linear  algebra  codes  and  sparse 
matrix,  algorithms  for  the  CRAY-1  (AFOSR/AFFDL) A  study  of 
the  relatively  poor  performance  of  a  sparse  matrix  code 
developed  without  the  aid  of  simulation  [8]  indicated  the 
importance  of  data  flow  considerations  in  addition  to  vector- 
izatxon  in  the  development  of  linear  algebra  algorithms  for  a 
memory-hierarchial ,  f unctionally-concurrent  processer  of  the 
CRAY-1  class.  Equation- solving  kernels  (full  solvers,  band 
solvers,  tridiagonal  solvers)  were  developed  and  timing  models 
produced  using  simulation  ,( 2 ][ 6 J .  Speedups  of  1.5:1  to  2:1 
viz-a-viz  LASL-prepared  assembly  language  codes  were  obtained. 
These  kernels  will  be  made  into  an  equation-solving  library 
for  the  CRAY-1  for  use  in  the  public  domain.  They  will  also 
be  incorporated  in  an  improved  sparse  matrix  algorithm  which 
has  asymptotic  speeds  approximately  four  times  that  of  the 
previously  developed  code. 

II.  Coupling  Activities 

A.  Air  Force/Government.  Besides  three  visits  to  AFFDL  to 


coordinate  the  aerodynamic  code  development,  seminars  were 
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presentea  to  Argonne  National  Laboratory  (June,  1978) 
and  to  Los  Alamos  Scientific  Laboratory  (May,  1973)  on  the 
topics  of  the  modeling  and  use  of  vector  processors. 

B.  Industrial .  Consultive  studies  with  Electric  Power  Research 
Institute  through  a  contract  with  General  Electric  Company 
were  completed.  These  involved  study  of  the  application 

of  vectorized  sparse  matrix  codes  to  electric  power  system 
analysis  and  produced  a  report  [7).  A  study  was  completed 
with  Mobil  Oil  Corporation  to  vectorize  computation  kernels 
used  to  solve  3-D  diffusion  equations  associated  with  oil 
reservoir  drilling  and  management. 

C.  Short  Course.  An  annual  one-week  short  course  entitled  High 
Speed  Computation:  Vector  Processing  was  continued  in  1978 
to  survey  the  state-of-the-art  in  vector  architectures  and 
related  problem-solving  metnods.  Among  the  attendees  were 
representatives  from  Rome  Air  Development  Center,  Air  Force 
Flight  Dynamics  Laboratory,  Air  Force  Weapons  Laboratory, 
and  Air  Force  Headquarters,  Pentagon. 
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